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ABSTRACT 


Sets  of  classical  turning  points,  together  with  the 

corresponding  vibrational  energies,  are  given  at  an  interval 

of  one  half  a  unit  of  quantum  number  for  the  following  molecular 

states:  C., ,  ,  k\  ,  b\  ,  c'H  ;  0_,  X^L",  B^f” ;  OH,  , 

^  u  g’  g’  g’  2’  g’  u’  i’ 

SiO,  XT"^,  A'n. 

These  numbers,  representing  points  on  potential  cunres, 
were  derived  from  Klein's  formulae  by  analytical  differentiation, 
and  numerical  integration  with  a  necessary  analytic  approximation 
over  a  small  part  of  the  range. 


1.  imROSUCTION 


An  earlier  Scientific  Report  (Jarmain  1960a)  contained 
tables  of  classical  turning  points  for  fifteen  diatomic  molecular  states, 
computed  by  a  modified  Klein-Dunham  series  method  described  in  detail 
elsewhere  (Jarmain  1960b).  At  higher  levels  of  unfavourable  states, 
however,  these  series  may  diverge,  or  converge  so  slowly  that  there  are 
not  enough  terms  formulated  to  produce  sufficiently  accurate  results. 

It  was  not  considered  worthwhile  to  work  out  coefficients  for  further 
terms  because  they  become  algebraically  unwieldy. 

The  procedure  adopted  was  to  perform  necessary  different¬ 
iations  einalytically,  thf'n  to  evaluate  integrals  numerically  on  an  I.B.M. 
650.  This  is  an  improved  method  of  constructing  potential  energy 
functions  up  to  a  vibrational  quantum  number  of  thirty.  The  only 
important  limitation  remaining  is  that  spectroscopic  obse3m/ations  are 
not  available  for  more  than  the  lowest  few  levels  of  certain  molecular 
states. 


i-CTHOD 


Tor  the  minim.um.  and  m.aximumi  values,  r.  „  of  the  internuclear 

*  t‘- 

distance  for  a  molecule  vibrating  with  an  energy  U,  Klein  (l932)  gives 
the  formula 


\j±th  f  and  g  defined  in  terms  of  the  auxiliary  function 


(2)  S(U,k)  = - ^  /[U-E(I,K)]  dl 

Ti(2u)  ' 


Here,  E(I,k)  is  the  sum  of  vibrational  and  rotational  energies  for  any 
level  up  to  U;  the  action  variable  I  =  h(v  +  l/2){  K  =  [j(j+l)h^  ]/871^H 
with  J  the  rotational  quantum  number;  p  is  the  reduced  mass  of  the 
molecule;  1=1'  when  E  =  U. 


is 

yii]K=o 


Replacing  v  +  l/2  by  V  and  dif ferentiating  under  the  integral  sign, 

-1/2 

' I  &  -  " 

U)  ° 

v 

0 

evaluated  at  K  =  0  for  the  normal  potential  cur^/e  (no  rotation). 

As  in  previous  work,  the  vibrational  quantum  number  is  not 
formally  restricted  to  integral  values.  The  integrals  in  (4)  cannot 
however  be  calculated  over  the  -whole  range  0  to  V  by  Simpson’ s  or 
other  method  of  numerical  integration,  because  the  integrands  become 
infinite  at  the  upper  limit  V .  This  difficulty  is  circumvented  by 
using  a  good  approximation  to  ^  near  V  =  V  ,  as  outlined 


The  program  written  for  the  I.B.M.  65O  permits  subdivision 


below . 
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of  one  unit  in  V  for  numerical  integration  according  to  tho  relation 

whore  H  is  tho  size  of  a  subdivision,  that  is  —  is  the  number  of  sub- 
divisions  per  unit  quantum  number,  is  the  maximum  V  to  bo  used,  and 
the  limit  of  519  is  dictated  by  computer  storage  space.  Suppose  for 
illustration  that  1 /H  is  chosen  to  be  24  S-na  V  =  3.5*  Energies  can 
be  interpolated,  if  not  already  known,  at  the  85  points  of  division 
using  arcs  of  two  parabolae.  Then  Simpson's  Rule  can  bo  applied  to 
tho  first  82  subdivisions,  leaving  tho  f-integrand  to  be  approximated 


in  tho  remaining  two  subdivisions  by  an  analytically  integrable  function. 

Tho  function  chosen  is  of  the  form  y  =  ax  ^  vfith  x  =  0  at 
V  =  V ,  _and  passing  through  (K,y^)  and  (2H,y2),  y^  and  y^  being  the 
second  and  third  last  ordinates  respectively.  Those  conditions  determine 


'a'  and  p  (tho  latter  is  very  near  one  half).  The  'end  contribution* 
to  tho  integral  is  easily  found  to  bo  ^Hy^/l-p.  with  1-p  approximated 
closely  by  tho  polynomial  2.6640  -  2.0403R  +  0.36067  R_  (from  a  Taylor 


Series  expansion  of  a  logarit’am) ,  whore  R  =  7  ^1^2' 

Includod  also  in  tho  integrand  for  g  (U)  is  -7^  ^  which 

d  nj  H 

«  2  2 

equals  q  B,  (Jarmain  1960b)  whore  q  =  u/h  ,  a  constant  that  can  be 

V  ^  ” 


combined  with  that  in  front  of  the  integral.  Three  B  values  are 
known  by  interpolation  for  the  'end  contribution'  to  g,  but  there  is 

h- 

t  % 

tho  question  of  how  they  should  be  weighted  when  tho  area  is  found 

2 

analytically.  If  a  parabola  y  =  Ax  +  Bx  +  C  is  passed  through  those 
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final  three  'points'  of  a  set  of  its  coefficients  can  be  evaluated 
and  the  produ  t  of  the  parabolic  function  with  y  =  ax  ^  integrated. 

The  weighting  factors  turn  out  to  be  1/15»  8/l5  s-nd  6/15  respectively, 
beginning  with  the  third  last  B^,  if  wo  take  p  to  be  one  half  for  this 
purpose.  A  weighted  B_^  is  calculated  accordingly  and  multiplied  into 
the  f-'end  contribution'. 

For  U,  B  and  B^^  in  wave  numbers,  the  appropriate  constant  in 

1/2  -1-1/2 

front  of  the  integral  for  f  is  r  B  ■  and  forgisr  B 
”  e  e  e  c 

These  constants  are  readily  included  in  the  program  so  that  f  and  g 
are  produced,  then  r^  ^  frof''  equation  (1). 

3 .  ^BuiS^d^To 

Turning  points  are  shown  to  three  decimal  places  at  an 

interval  of  one  half  in  V  in  tables  "'-10  for  the  follovring  molecular 

states;  C,,  ,  A^n  ,  3^0  ,  c'H  j  0.,  B^Z";  OH,  X^n.  , 

2’  _  u’  g’  g  g’  2’  g’  u’  ’  1’ 

GiO,  x'^Z'*',  A. "'ll.  Potentials  for  the  ground  states  of  C2  and  OH  were 
previously  published  (jarmain  1960a)  but  have  been  extended  and  revised 

slightly.  It  is  perh^s  worth  mentioning  that  agreement  is  nearly 

^  £ 

perfect  to  four  figures  between  I.B.X.  results  for  certain  states  not 

listed  in  this  report,  and  those  published  last  year  on  the  basis  of 
*  m 

the  series  method. 

Bxperimental  data  have  been  more  carefully  surveyed  than 
before,  reference  having  been  made  in  almost  every  instance  to  original 
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papers  named  by  Herzberg  (195'-^)  ‘'ind  papers  published  since  1950*  The 
source  of  the  data  is  given  with  each  table. 

Table  6  deserves  special  comment  since  turning  points  have 
been  calculated  to  the  level  v  =  21,  which  is  within  possibly  12-13  cm  ^ 
of  the  convergence  limit  (Brix  and  Herzberg  1954).  i^valucs  on  the 
repulsive  side  of  the  curve  reach  a  minimum  of  1.330A  at  about  V  =  l7.5i 
then  show  an  unreal  increasing  trend,  varying  irregularly  from  1.33l 
to  1.338.  Such  behaviour  is  likely  due  to  high  sensitivity  of  f  and  g 
to  small  errors  in  the  experimental  data  in  that  region,  because  these 
produce  large  changes  in  the  function  ;3(cf.  Vanderslice,  Mason  and 
M.aisch  I960)  • 

APPLICATION 

The  intention  is  to  use  tabulations  of  'true'  potential 
energy  functions  U(r),  computed  in  this  way  (but  inverted  by  interpolation 
so  that  r  is  the  independent  variable) ,  in  the  reduced  Schrftdinger 
equation 

(5).  ^  ^  [b-  U(r)  ]  .  .  0  , 

dr^  h.  ^ 

Methods  of  solving  (5)  Tor  more  accurate  vibrational  wave  functions 
by  step-by-step  integration  processes  arc  at  present  being  tested. 


TABLE  1 


2  u 


V 

Energy 
(cm  ') 

I’l  (a) 

r2(A) 

O.B 

817*8 

1.257 

1 .374 

1 .0 

1629.7 

1.236 

1.402 

1-5 

2435.8 

1.220 

1.425 

2.0 

3236.1 

1.208 

1  •  445 

2.5 

4030.4 

1.197 

1.463 

3.0 

4819. 0 

1.188 

1.479 

3.5 

5601.8 

I.I79 

1-495 

4*0 

6378 .7 

1 . 1 7 1 

1.510 

4.5 

7149.9 

I.I64 

1.525 

5.0 

7915.0 

1.15s 

1.539 

5-5 

.  8674.3 

1.152 

1-553 

6.0 

9427.9 

1 .  146 

1 .566 

6.5- 

IOI76 

1 . 140 

1.579 

7.0 

109 18 

1.135 

1.592 

7-5 

11654 

1.130 

1 .605 

8.0  • 

12384 

1.126 

1.61.8 

3-5 

13108 

1.121 

1.630 

9.0 

1382  7 

I.II7 

1.642 

9.5 

14540 

1.113 

.  1.655 

10.0 

15247 

.1.109 

1.667 

10.5 

1594.8 

1.105 

1 .678 

11.0 

16643 

1.102 

1.690 

11.5 

I7332 

1.098 

Si  1.702 

"12  .Q, 

18016  • 

^  1 .095 

1.714 

J2.5 

I8693 

1.091 

*  1.726 

Phillips  (1943) 
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TABLE  2 
C„  A^n 


V  'Energy 


0.5 

889  *9 

1.0 

1771,5 

1-5 

2643«6 

2.0 

3506. 5 

2.5 

4360.0 

3.0 

5203.2 

3.5 

6035.6 

4.0 

6857.2 

4.5 

7667.9 

5.0 

8469.3 

5.5 

9253.6 

6.0 

10021 

6.5 

IO77I 

7  .0 

11507 

7-5 

12227 

8.0 

12929 

8.5  ■ 

13614 

9^0 

-  14277 

9.5' 

14919 

10.0 

15539 

10.5  ’ 

16138 

9 

Phillips  (1945) 


t-iCa) 

r2(A) 

1  *2  "i  3 

1 .325 

1.193 

1.352 

1.178 

1.375 

1 » 166 

1.394 

'1.156 

1  .412 

1.147 

1.429 

1.139 

1.445  • 

1.132 

1.461 

1.126 

1 .476 

1.120 

1.492 

1.114 

1.507 

1 . 109 

1 .523 

1.104 

1.539 

1 , 100 

■'•554 

1 . 096 

1.570 

1.092 

1-587 

1.088 

1 .603 

1.0P,4 

1 .620 

1.081 

1.638 

I.O77 

1.657 

_ w 

• 

0 

• 

1.676 

9 


TABLE  3 


B^n 


B 


V 

Energy 

_ _ 

r.,(A) 

0.5 

543.8 

1.470 

1 .614 

1.0 

IO7O.6 

1.445  ' 

1.652 

1-5 

158O.4 

CM 

• 

1.685 

2.0  • 

2075.9 

1.413 

1 .714 

2.5 

2559.9 

1.401 

1 .741 

3.0 

3032.4 

1.391 

1 .766 

3.5 

3493.3 

1.382 

1.791 

4.0 

3945.0 

1.373 

I.8I5 

4.5 

4388.2 

1.366 

1.838 

5.0 

4822.9 

1.359 

1.861 

5-5 

5249-0 

1.352 

1.884 

Phillips  (^949) 
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TAB13  4 

C.  c^n 

■=:  g 


V 

Energy 
(cm  ”') 

r/A) 

r2(A) 

0.5 

900.1 

1 .203 

1.315 

1.0 

1790. 8 

1.183 

1.342 

1-5 

2664.5 

1 .169 

1.365 

2.0 

3521,3 

1.157 

1.386 

2.5 

4361.1 

1.147 

1.405 

3.0 

5I8O.4 

1.139 

1.425 

3.5 

5965.7 

1 .131 

1.445 

4.0 

67I6  .9 

1.125 

1 .466 

4.5 

7433.9 

1.119 

1.488 

5.0 

8O96 .2 

1.114 

1.515 

5-5 

8693.0 

1.110 

1.545 

6.0 

9224 . 1 

1.106 

1.580 

6.5 

9689-7 

1.102 

1.620 

Herzberg  and  Sutton  (I94O) 
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5 

0^  X^Z" 
2  g 


V 

Energy 

(cm 

r^(A) 

t^{a) 

0.'5 

787.2 

1.159 

1 .262 

1 .0 

,  1568.4 

1 . 140 

1.287 

1-5 

2343.6 

1.127 

1.308 

2.0 

3112.8 

1.116 

1.325 

2.5 

3876.2 

1.106 

1.341 

3.0 

4633.8 

1.098 

1.356 

3.5 

-  5385-5 

1 .091 

1.370 

4.0 

6 131. 4 

1.084 

1 .384 

4.5 

667I.5 

1.078 

1.397 

5.0 

7605.9 

1 .072 

1 .410 

5-5 

8334 .6 

1 .067 

1.422 

6.0 

905^.5 

1 .062 

1.434 

6.5 

9774.- 

1.057 

1.446 

7  *0 

1048'^ 

1.053 

1.457 

'^•5 

’II92 

1.049 

1.469 

8.0 

1 1392 

1 .045 

1.480 

8.5 

12537 

1 .041 

1.491 

9.0 

'3^76 

1 .03? 

1 .502 

9.5 

13959 

1.034 

1.513 

10.0 

14637 

^  1.031 

1.524 

10.5 

15309 

1 .027 

1.534 

11.0 

'59-5 

-  1.024  ' 

1.545 

11.5 

16636 

1 .021 

1.556 

12.0 

I729O 

1 .013 

1 .566 

12.5 

..  17940 

1 .016 

1 .577 

13.0 

18583 

1 .013 

1.587 

13.5 

19221 

1 .011 

1."'8 

14.0 

19854 

1 .008 

1.608 

14.5 

20480 

1 .005 

1.619 

15.0 

21101 

1 .003 

1 .629 

15.5 

2I7I6 

1 .001 

1.640 

Feast  ( 1950) 

Herzberg 

(1950) 
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TABLE  6 
0,, 

2  u 


V 

Energy 
(cm  ^) 

ri(A) 

0.5 

343.1 

1.531 

1.688 

1 .0 

686.8 

1.505 

1.727 

1-5 

1031.1 

1  .'I87 

1.758 

2.0 

1366.7 

1 .47I 

1.788 

2.5 

1696.3 

1.458 

I.8I5 

3.0 

2020.0 

1.446 

1.840 

3.5 

2337.7 

1 .436 

1.865 

4.0 

2649.7 

1 .427 

1.889 

4.5 

2955.3 

1.419 

1 .912 

5.0 

3254.4 

1.411 

1.936 

5-5 

3547.1 

1.405 

1 .960 

6.0 

38  31.1 

1.398 

1.985 

6.5 

4103.3 

1.392 

2.008 

7.0 

4378.7 

1.386 

2.032 

7-5 

4642 .3 

1.380 

2.056 

fl.O 

4897.9 

1.376 

2.083 

5.5 

5144.1 

1.373 

2.111 

9.0 

53p,i.7 

1 . 369 

2.139 

9.5 

5609.9 

1.365 

2.168 

10.0 

5829.1 

1.360 

2 . 196 

10.5 

6  0  38  .8 

1.355 

2.226 

11.0 

62  38  .0 

1.351 

2.259 

11.5 

6427.0 

1.349 

2.296 

12.0 

6603.3 

'  1.346 

2.336 

12.5 

6770.1 

1.344 

2.377 

13.0 

6925.8 

1.341 

2.420 

13.5 

7O7I .0 

1.339 

2.467 

14.0 

7264.1 

1.337 

2.520 

14.5 

7326.0 

1.336 

2.577 

15.0 

7436.6 

1.334 

2.639 

15.5 

7  536  .0 

1.333 

2.709 

16.0 

7625.0 

1.332 

2.783 

16.5 

77O5.O 

1.331 

2.864 

I7  .0 

■7776.0 

1.331 

2.954 

17.5 

7337.9 

1.330 

3.057 

18.0 

7892 . 3 

3. 165 

18 .5 

7940.0 

3.285 

19.0 

7981.2 

3.420 

19.5 

3015.7 

Average 

3.576 

20,0 

8044.6 

1.335 

3.750 

20.5 

8068 . 3 

3.957 

21  .0 

8086.9 

4. 217 

21.5 

8100.3 

4.574 

Loch.te-Holtgreven  and  Dieke  (1929) 
Curry  and  Herzberg  ( 1934) 

Knauss  and  Ballard  (1935) 

Brix  and  Herzberg  (1954) 
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TABIE  7 

OH  ^n. 
1 


V 

Energy 

(cm 

r^(A) 

^2^  A) 

0.5 

1848  -5 

0.883 

1.080 

1.0 

3654.5 

0.852 

•1  4  *5  X 

1  •  1  jq. 

1-5 

54I8.I 

0.831 

1.179 

2.0 

7140.6 

O.8I4 

1 .219 

2.5 

8822 . 1 

0.799 

1.257 

3.0 

10463 

0.787 

1.294 

3.5 

12062 

0.776 

1.330 

4.0 

13621 

0.767 

1.365 

4.5 

I5I4O 

0.758 

•  1.399 

5.0 

16618 

0.751 

1.434 

5-5 

IRO55 

0.744 

1 .469 

6.0 

I945I 

0.737 

1.503 

6.5 

208 O7 

0.731 

1.538 

7.0 

22121 

0.726 

1.573 

7-5 

23395 

0.720 

1.609 

8.0 

24623 

0.715 

1 .646 

8.5 

258 O7 

0.7II 

1.684 

9.0 

26947 

O.7O7 

1.724 

9.5 

28042 

O.7O3 

1 .764 

Herman  and  Hombeck  (''95^) 
Chamberlain  and  Roesler  (1955) 
Wallace  ( I96O) 
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TABLE  8 


V 

Energy 

(cm 

r^(A) 

0.5 

1566.5 

O.9I8 

1.133 

1.0 

3084.8 

0.886 

1.193 

1-5 

4555.1 

0.863 

1 .244 

2.0 

5976.5 

0.846 

1.292 

2.5 

7348.0 

0.832 

1.338 

3.0 

8669.6 

O.8I9 

1.383 

3.5 

9941.3 

0.809 

1.427 

4.0 

11162 

0.800 

I.47I 

-5 

12331 

0.791 

1 .5I6 

5.0 

13448 

0.784 

1.562 

5*5 

I45I4 

0.777 

1.609 

Bieke  and  Crosswhite  (1948) 
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TABLS  9 
SiO 


V 

Energy 

(cm 

r^(A) 

1-2  (a) 

0.5 

619.2 

1 .460 

1.563 

1.0 

1235.5 

1.441 

1.588 

1.5 

1848-8 

1 .427 

1 .607 

2.0 

2459.2 

1  .416 

1 .624 

2.5 

3066.6 

1 .406 

1.639 

3.0 

3671. 0 

1.397 

1.653 

3.5 

4272.5 

1.389 

1 .666 

/.o 

487 1.0 

1 .38,2 

1.679 

4.5 

5466 .6 

1.375 

1.691 

9.0 

6059  -2 

1.369 

1.703 

5-5 

6648-8 

1 , 364 

1.714 

6.0 

7235.5 

1.358 

1.725 

6.5 

78I9.2 

1.353 

1.736 

7.0 

8400.0 

1 . 348 

1.746 

7-5 

8977-8 

1 . 344 

1.756 

8.0 

9552.6 

1.339 

1.766 

8.5 

10125 

1.335 

r.776 

La^erq-/:  -t 

and  Uhler  (1953) 
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TABLE  10 

SiO  aVi 

V 

Energy 
(cm  ^ 

r^(A) 

t^{a) 

0.5 

424.7 

1.562 

1.637 

1.0 

846.3 

1.540 

I.7I7 

1 .  c; 

-  J 

•iOCf.  c 

•  t.  i .  %  \J 

1.524 

1-742 

2.0 

1679.7 

1.511 

1.763 

2.5 

2091.5 

1 .500 

1.783 

3.0 

2500.2 

U49O 

1.801 

3.5 

2905.6 

1.481 

1.8 19 

Lagerqvist  and  Uhler  (1953) 
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